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ABSTRACT 

Continuous  Phase  Frequency  Shift  Keying  (CPFSK)  is  a 
potentially  attractive  modulation  scheme  with  constant 
envelope  and  good  spectral  characteristics  for  energy 
constrained  and  band- limited  satellite  channels. 

This  research  deals  with  Orthogonal  Frequency  Shift 
Keying  (FSK) ,  Minimum  Shift  Keying  (MSK)  which  is  a  special 
case  of  CPFSK,  uncoded  quaternary  CPFSK  and  finally  coded 
quaternary  CPFSK.  Orthogonal  FSK  is  simulated  by  making  the 
modulation  index  (h)  equal  to  one,  and  all  the  other 
simulations  are  performed  with  h=l/2 . 

A  rate  1/2  convolutional  encoder  with  constraint  lengths 
(k)  ,  k=2,  3,  and  4  are  used  in  coded  quaternary  CPFSK 
simulations.  Good  coding  gains  are  obtained  with  only  a  slight 
increase  in  receiver  complexity. 

Soft  decision  with  the  Viterbi  Algorithm  was  applied  to 
all  CPFSK  and  one  MSK  application.  Hard  decision  was  applied 
to  Orthogonal  FSK  and  another  MSK  application. 
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I .   INTRODUCTION 

An  increasing  need  for  communication  between  individuals 
and  nations  is  making  satellite  communication  more  and  more 
attractive  every  day.  Physical  constraints  of  space  based 
applications  force  us  to  choose  bandwidth  and  energy 
efficient  modulation  techniques. 

Continuous  Phase  Modulation  (CPM)  is  a  digital  modulation 
scheme  with  continuous  phase  and  constant  envelope.  Continuous 
phase  makes  it  spectrally  efficient,  and  constant  envelope 
seems  desirable  where  non-linearities  are  present  in  space 
segment  power  amplifiers  and  repeaters. 

Coding  is  used  to  reduce  reception  errors.  The  trade-off 
in  coding  is  a  decrease  in  bandwidth  efficiency.  This  thesis 
mainly  presents  the  research  for  rate  1/2  convolutional 
quaternary  CPM  schemes  with  different  constraint  lengths. 

The  starting  point  for  this  research  was  the  simulation  of 
rate  2/3  trellis  coded  CPFSK.  The  complexity  and  apparent  poor 
performance  of  this  code  forced  us  to  first  study  orthogonal 
FSK,  and  then  continue  with  MSK,  uncoded  quaternary  CPFSK  and 
finally  finish  with  rate  1/2  trellis  coded  CPFSK.  Due  to  the 
limitation  of  time,  we  were  unable  to  return  to  the  2/3  rate 
codes . 

The  simulations  were  done  using  PRO-MATLAB  software  run  on 


Sun  SPARC  station  computers.  Hardware  limitations  and  the 
complexity  of  the  CPM  schemes  made  this  research  difficult  and 
very  time  consuming. 

Chapter  II  gives  a  brief  review  of  Continuous  Phase 
Modulation.  Chapter  III  focuses  on  Orthogonal  FSK,  and  two 
different  receiver  structures  for  MSK.  Chapter  IV  presents 
uncoded  quaternary  CPFSK  and  coded  quaternary  CPFSK  with 
different  constraint  lengths.  Finally,  Chapter  V  focuses  on 
the  final  conclusions  and  opportunities  for  future  research. 

An  identical  string  of  fifty  thousand  information  bits  has 
been  used  for  each  simulation  to  provide  a  comparison  between 
systems . 


II.  BACKGROUND 

This  chapter  focuses  on  the  basic  characteristics  of 
Continuous  Phase  Frequency  Shift  Keying.  In  CPFSK,  frequency 
shifts  are  used  to  transmit  information.  The  phase  of  the 
CPFSK  signal  is  continuous  during  frequency  shifts.  This 
prevents  the  large  spectral  side  lobe  effects  that  are 
observed  in  FSK,  and  PSK.  There  is  also  memory  in  CPFSK  due  to 
the  continuous  phase.  Phase  trellises  will  be  presented  in 
this  chapter  later  on. 

A.   CPFSK  SIGNAL  DESCRIPTION 

Equation  (2.1)  represents  the  transmitted  signal  in  CPFSK, 
[Ref.  1] 


s(t)  = 


2^cos[2*fct+0(t;I)+0o]  (2'1) 


where  Ts    is  the  symbol  duration,  Es    the  energy  per  symbol, 

9(t;J)  is  the  time  varying  phase  modulation  induced  by  the 
input  data  sequence,  fc   is  the  carrier  frequency,  and  finally,  0O 

is  the  initial  phase  of  the  carrier  assumed  to  be  zero  for 
simplicity.  Equation  (2.2)  represents  the  phase  of  the  carrier 
in  the  interval  nT„<  t<  (n  +  1)  Tc . 


e(t;  I)  =*££  Ik+2%hq{  t-nTs)  (2.2) 


The  parameter  h   is  called  the  modulation  index.  Typically  h   is 

chosen  to  be  a  rational  number  h  =  —     where  in  and  p  are 

P 


relatively  prime  positive  numbers,  in  order  that  9(t; J)  takes 
on  a  finite  number  of  states.  The  q(t)  is  a  phase  shaping 
pulse  signal,  I  is  the  input  signal.  For  M-ary  CPFSK,  M=2n 
and  In    is  one  of  the  {±1,  ±3, ,+(M-l)}  information 

symbols . 

Phase  shaping  signals  can  be  in  a  rectangular  or  raised 
cosine  form.  A  rectangular  form  is  used  in  this  research. 
Phase  trajectories  for  rectangular  pulse  shapes  are  piecewise 
linear.  There  are  p  phase  states  when  m  is  even  and  2p  phase 
states  when  in  is  odd.  Smoother  phase  trajectories  can  be 
obtained  by  using  raised  cosine  pulse  shapes. Phase 
trajectories  for  binary  CPFSK  are  illustrated  in  Fig.  2.1. 

B.   ERROR  PERFORMANCE 

Minimum  Shift  Keying  (MSK)  was  used  as  a  baseline 
modulation  scheme;  a  special  case  of  binary  CPFSK  with  h=l/2. 
In  coded  CPFSK  cases  we  also  made  comparisons  with  MSK, 
because,  although  we  send  two  bits  at  a  time,  one  of  them  is 
an  information  bit,  and  the  other  is  a  redundant  coding  bit. 

Equation  (2.3)  represents  the  probability  of  error  in  the 


Fig.  2.1:  Phase  trellis  for  binary  CPFSK, 


presence  of  additive  white  Gaussian  noise  for  MSK.  [Ref.  2] 

(2.3) 


Pe=0{, 


2^ 
N  N0 


The  error  performance  can  be  improved  by  increasing  the 
signal  energy,  but  in  satellite  applications  there  are 
limitations  on  power  sources. 

Error  performances  for  each  modulation  scheme  are  given  at 
the  end  of  each  section. 


C .   BANDWIDTH 

Bandwidth  is  the  range  of  frequencies  which  contains  99 
percent  of  the  total  signal  power.  A  smaller  value  of  h  will 
increase  the  bandwidth  efficiency  with  an  increase  in 
probability  of  error. 


III.  FSK-MSK 

This  chapter  focuses  on  orthogonal  FSK  and  MSK.  The  same 
modulator  and  bit  mapper  was  used  in  all  simulations.  Input 
bits  were  mapped  to  modulation  symbols  according  to  Table  3.1. 

TABLE  3.1:  MAPPING  FOR  ORTHOGONAL  FSK  AND  MSK. 

input  bit      mapper  (I) 

0  >  -i 

1  >   1 

A.   ORTHOGONAL  FSK 

Orthogonal  FSK  simulation  is  accomplished  by  using  binary 
CPFSK  with  h-1.  Coherent  detection  is  used  in  the  demodulator 
and  hard  decisions  are  made  to  detect  the  received  signal.  By 
making  h=l ,  we  have  two  phase  states;  0  and  rc.  The  phase 
trellis  for  this  system  is  illustrated  in  Fig  3.1.  The  odd 
symbols  represent  the  system  at  state  0,  and  even  symbols 
represent  the  system  at  state  n  .  Multiplying  the  received 

signal  by  ejnt/Ts      offsets  the  total  phase  to  0  at  the  end  of 

each  symbol  and  decreases  the  receiver  complexity.  The 
maximum  value  of  the  correlator  outputs  for  the  two  modulation 
symbols  is  selected  as  the  decision  variable. 


Fig.  3.1:  Phase  trellis  for  Orthogonal  FSK 


The  MATLAB  source  code  used  for  this  system  is  listed  in 
Appendix  A.  Equation  (3.1)  represents  the  probability  of  bit 
error  for  orthogonal  FSK.  [Ref.  2] 


Pe=Q(. 


(3.1) 


N  ^o 


The  performance  of  this  system  is  illustrated  in  Fig.  3.2.  The 
orthogonal  FSK  simulation  was  done  for  simulation 
verification,  and  as  guidance  to  MSK  simulation. 


Fig.  3.2:  Performance  of  Orthogonal  FSK 


B.   MSK  WITH  VITERBI  ALGORITHM  DECODING 

MSK  is  a  special  case  of  binary  CPFSK  with  h=l/2  where  the 
phase  changes  ±n/2  at  every  symbol  (see  Fig.  3.3)  .  The  total 
phase  is  +rc/2  for  odd  symbols  and  the  total  phase  is  0  or  n 
for  even  symbols.  When  the  received  signal  is  multiplied  by 
ej2*t/4r%  it  offsets  the  phase  such  that  the  total  phase 
changes  either  by  0  or  by  it  at  every  symbol .  This  reduces  the 
number  of  states  required  in  the  Viterbi  decoder.  [Ref.  3] 


Fig.  3.3:  Phase  Trellis  for  MSK 

The  state  trellis  for  the  MSK  waveform  is  illustrated  in 
Fig.  3.4.  There  are  two  states,  0  and  tc  .  There  are  two  paths 
out  of  each  node  corresponding  to  the  two  possible  values  of 
the  input.  Information  bit  1  produces  a  change  in  state  and 
information  bit  0  produces  no  change.  The  Euclidean  distances 
of  the  transmit  symbol  waveforms  from  the  received  signal  are 
used  as  input  to  the  Viterbi  algorithm;  the  Viterbi  algorithm 
chooses  the  path  with  the  minimum  total  distance. 
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State 
1 


0 


Fig.  3.4:  State  Trellis  for  MSK 

The  transition  matrix  of  the  MSK  with  Viterbi  algorithm 
decoding  is  illustrated  in  Table  3.2. 


TABLE  3.2:  VITERBI  TRANSITION  MATRIX. 


10  2  2  13 

1112  0  4 


Each  row  of  the  matrix  corresponds  to  a  state.  The  first 
column  of  the  matrix  represent  the  states  that  we  are  coming 
from,  the  second  column  represents  the  input  value,  and  the 
third  column  represents  the  corresponding  number  of  the  point 
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in  the  signal  constellation.  Columns  4,  5,  6  are  the  same 
respectively.  The  MATLAB  source  code  for  this  system  is  given 
in  Appendix  B . 

Theory  vs.  simulation  of  the  MSK  bit  error  probability 
with  Viterbi  detection  is  illustrated  in  Fig.  3.5. 
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Fig.  3.5:  Performance  of  MSK  with  Viterbi  Detection. 

C.   MSK  WITH  DETECTION  USING  TWO  SYMBOLS  AT  A  TIME 

Equation  (2.1)  can  be  rewritten  as  (3.2).  [Ref.  4; 
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)  =  2^cos[e(t)  Jcos(2afct)  - 

M  Tb      : (3.2) 

\  2— 2 sin  [  0(  t)  ]  sin(27ifct) 


Equation  (3.3)  shows  that  the  phase  of  an  MSK  signal  increases 
or  decreases  linearly  within  the  bit  duration.  [Ref.  4] 

0(t)  =0(o)Tu-|-  (3.3) 

2Tb 

The  minus  sign  in  (3.3)  represents  information  signal  0,  and 
the  plus  sign  represents  one.  The  signal  phase  can  have  values 
of  ±tc/2  at  odd  multiple  of  Tb   or  0  and  tc  at  even  multiple  of 

Tb    (see  (3.3)  and  Fig.  3.3). 

When  (3.3)  is  inserted  in  (3.2),  it  is  clear  that  the 
input  sequence  during   (0,7^)   affects  only  the  quadrature 

component  of  (3.2),  and  not  the  in-phase  component  over  the 
signalling  interval   (0,2  7^)  .  Similarly  the  input  sequence 

during  {Tb,2Tb)    affects  only  the  in-phase  component  of  (3.2), 

and  not  the  quadrature  component  over  the  signalling  interval 
(T,3Tb)  . 

Table  3.3  was  constructed  from  the  information  discussed 
above.  [Ref.  4] 


13 


TABLE  3.3:  EFFECT  OF  INFORMATION  BITS  ON  PHASE 


INFORMATION   BITS 

9(0) 

®(Tb) 

0 

71 

tc/2 

0 

0 

-n/2 

1 

0 

7t/2 

1 

7t 

-It/2 

The  hard  decision  in  the  demodulator  is  done  according  to 
Table  3.3,  and  the  MATLAB  source  code  for  this  system  is 
illustrated  in  Appendix  C. 

Theory  vs.  simulation  of  MSK  bit  error  probability  for 
detection  using  two  symbols  at  a  time  is  illustrated  in  Fig. 
3.6. 


Fig.  3.6:  Performance  of  MSK  for  detection  using  two 
symbols  at  a  time. 
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IV.  CODED  CPFSK  SCHEMES 

Channel  coding  can  be  used  to  reduce  the  number  of  errors 
when  a  message  is  transmitted  over  a  noisy  channel.  To  achieve 
this  goal,  redundancy  is  added  to  the  information  sequence. 

The  encoder  used  in  this  research  is  a  linear,  sequential, 
finite  state  machine  which  is  implemented  by  using  shift 
registers.  Rate  1/2  encoders  with  constraint  lengths  of  2,  3, 
and  4  are  used.  The  transmitter  block  diagram  for  these 
systems  is  illustrated  in  Fig  4.1,  and  receiver  block  diagram 
is  illustrated  in  Fig.  4.2. 


SOURCE 


CONVOLUTIONAL 
ENCODER 
RATE  1/2 


MAPPER 


M-ary  CPM 
MODULATOR 


OUT 

3 


Fig.  4.1:  Transmitter  Block  Diagram. 
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SOURCE 

DECODER 


OUT 

l=> 


Fig.  4.2:  Receiver  block  diagram. 

The  bit  mapping  used  in  this  chapter  is  illustrated  in 
Table  4.1.  Two  coded  bits  are  sent  in  each  symbol  duration. 


TABLE  4.1:  MAPPING  FOR  CODED  4-CPFSK. 

input  dibits      mapper  (I) 

00   >   -3 

0i   >   -1 

10   >    l 

11  -- --->    3 


The  phase  trellis  for  the  uncoded  quaternary  CPFSK  with 
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h=l/2  is  illustrated  in  Fig.  4.3.  By  combining  coding 
techniques  presented  in  this  chapter,  the  probability  of 
choosing  the  wrong  information  sequence  will  be  reduced 
significantly. 


Fig.  4.3:  Phase  trellis  for  quaternary  CPFSK, 


The  demodulator  used  in  this  chapter  is  illustrated  in  Fig 
4.4.  As  is  seen,  the  received  signal  is  first  offset 
and  then  passed  through  four  correlators.  Equation  (4.1) 
represents  the  correlator  outputs. 


Cn=jrfQTsz*(t)  .sn(t)dt 


(4.1) 
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Fig.  4.4:  Demodulator. 

Real  parts  of  the  correlator  outputs  are  taken,  and  are  used 
as  elements  of  a  four  dimensional  received  signal  vector.  The 
euclidean  distances  of  this  vector  are  calculated  in  the 
comparator  from  each  of  the  vectors  in  the  signal 
constellation.  The  vectors  of  the  signal  constellation  are 
determined  as  follows. 

In  the  absence  of  noise,  the  received  signal  will  be 

r(t)  =eij9k)  e{JKhI*t/Ts)    where  the  phase  state  &k=(0,K)     and  the 

transmitted  information  sequence  Ik=  ( -3  ,  -1, 1 ,  3  )  .  Thus  there 
are  eight  vectors  in  the  signal  constellation  determined  by; 
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Ck=  —  Re 


[f*sr*(t)sn(t)dt] 
Ckn  =  jrRe  [  e  lj»*]  f  T'etJ%htz^  t/T*dt] 


^kn     T kn 


(4.2a) 

(4.2b) 

(4.2c) 
(4. 2d) 


Table  4.2  illustrates  the  T kn    in  a  tabulated  form.  Each  row 

of   Table   4.2   represents  a  single  point   in  the  signal 
constellation. 

TABLE  4.2:  SIGNAL  CONSTELLATION. 


v^a 

-3 

-1 

1 

3 

-3     (0^=0) 

1 

0 

0 

0 

-3     (@k=n) 

-1 

0 

0 

0 

-1     (0^=0) 

0 

1 

0 

0 

-1       (0^=71) 

0 

-1 

0 

0 

1     (9^=0) 

0 

0 

1 

0 

1     I9k=it) 

0 

0 

-1 

0 

3     (0^=0) 

0 

0 

0 

1 

3     (ek=n) 

0 

0 

0 

-1 

Equation  (4.1)  is  represented  as  (4.3)  in  the  MATLAB 
program. 


TjAt 

Cn=  £  r*(kAt)   SjkAt) 


(4.3) 
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where  At  is  the  sampling  period,  and  rs/At  is  the  number  of 

samples  per  symbol   (N)  .   In  MATLAB  simulation,   Table  4.2 
becomes  Table  4.3. 

TABLE  4.3:  SIGNAL  CONSTELLATION  FOR  MATLAB  APPLICATION. 


V^n 

-3 

-1 

1 

3 

-3     (6^=0) 

N 

0 

1 

0 

-3     (@k=n) 

-N 

0 

-1 

0 

-1     (9^=0) 

0 

N 

0 

1 

-1       (0^  =  71) 

0 

-N 

0 

-1 

1     (Qk=0) 

1 

0 

N 

0 

l    (ek=n) 

-1 

0 

-N 

0 

3     (8^=0) 

0 

1 

0 

N 

3     {Qk=n) 

0 

-1 

0 

-N 

A.   CODED  4-CPFSK  WITH  CONSTRAINT  LENGTH  TWO 

The  coding  used  in  this  section  is  created  by  using  two 
shift-registers.  [Ref.  5]  The  encoder  is  illustrated  in  Fig. 
4.5,  and  the  state  trellis  for  this  system  is  illustrated  in 
Fig  4.6  (b) . 

The  demodulator  calculates  the  euclidean  distances  of  the 
received  signal  vector  from  the  eight  vectors  illustrated  in 
Table  4.4.  At  each  symbol  duration,  eight  distance  values  are 
entered  to  the  Viterbi  algorithm  and  the  true  path  is  chosen 
according  to  minimum  distance. 
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Most  significant  bit 


Least  significant  bit 


Fig.  4.5:  Encoder, 


States 


1 

-3 

^^ 

P\ 

2 

/*      \ 

(a) 


(b) 


Fig.  4.6:  Phase  Trellis  (a),  and  State  Trellis  for  coded  4 
CPFSK  with  k=2  (b) . 
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When  comparing  Fig.  4.3  and  Fig  4.6  (a) ,  it  is  clear  that 
at  each  symbol  duration,  the  probability  of  choosing  the  wrong 
state  is  reduced  by  this  simple  coding.  Later,  in  the 
following  sections,  more  complicated  codes  will  be  studied  and 
the  coding  gains  that  are  obtained  will  be  seen  more  clearly. 


TABLE  4.4:  SIGNAL  CONSTELLATION. 


0,-1,0,  -AT  0,1,  0,N     l,0,N,0       -1,0, -N,  0 
0,-^,0,-1  0,N,  0,1  N,  0,1,0  -N,  0,-1,0 


The  transition  matrix  for  this  system  is  illustrated  in  Table 
4.5. 


TABLE  4.5:  VITERBI  TRANSITION  MATRIX. 


101102213214 
125126237238 


Each  row  of  the  matrix  corresponds  to  a  state.  The  first 
column  of  the  matrix  represent  the  states  that  we  are  coming 
from,  the  second  column  represents  the  input  value,  and  the 
third  column  represents  the  corresponding  point  in  the  signal 
constellation  that  is  illustrated  in  Table  4.4;  the  points  are 
numbered  rowwise.  All  other  columns  are  the  same  respectively. 
The  performance  of  this  system  is  illustrated  in  Fig.  4.7. 
The  coding  gain  obtained  will  not  be  worth  what  must  be 
sacrificed  in  bandwidth  efficiency  (see  Fig.  4.7).  Since  this 
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was  the  initial  experiment  in  the  use  of  coding,  it  provided 
a  baseline  for  our  study.  The  MATLAB  source  code  for  this 
system  is  listed  in  Appendix  D. 


-OOOO.l.ty      for     MSK      &r>a      4-CPFSK      ~,tn      k  — 2 


Fig  4.7:  Performance  of  the  coded  4-CPFSK  with  k=2 . 

B.   CODED  4-CPFSK  WITH  CONSTRAINT  LENGTH  THREE. 

The  coding  used  in  this  section  was  created  by  using  three 
shift-registers.  The  encoder  is  illustrated  in  Fig.  4.8,  and 
the  state  trellis  for  this  system  is  illustrated  in  Fig  4.9. 
[Ref.  5] 

The  demodulator  calculates  the  euclidean  distances  of  the 
received  signal  vector  from  the  eight  vectors  presented  in 
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Most  significant  bit 


Least  significant  bit 


Fig   4.8:    Encoder. 


States 


00 

-3 

01 

.j^CS 

10 

?c?^ 

11 

^"^^ 

Fig  4.9:  State  Trellis  for  coded  4-CPFSK  with  k=3 
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Table  4.4.  At  each  symbol  duration,  eight  distance  values  are 
entered  to  the  Viterbi  decoder  and  the  true  path  is  chosen 
according  to  minimum  distance. 

The  transition  matrix  for  this  system  is  given  in  Table 
4.6. 
TABLE  4.6:  VITERBI  TRANSITION  MATRIX. 


101102225226 
337338413414 
202201126125 
314313438437 


The  performance  of  this  system  is  illustrated  in  Fig. 
4.10. 


tlTOr       =-"0000.1 

ty      'or      MSK      or,  a      4-CPFSK      with      K-     3 

; j 

-T-T-r-T        ._^~^-r^i_^__ 

^--^iiH^U^^ 

[ 

iillliiiH 

! !!!!!!!: !!!!!!!!!  11 !!!! ! 

\  :  Go^'^y  :  !i  :  1  :  :  :  :^K<|H 

-•Tr.W<}^--    ;  •  " 

||  { ! 1 ! I !  { i 

!::!!!!;!!:  1 !!!!  !i !! 1 !  : 

:  :                         !  ■      '         '     '        ; 

1!!!!!!!!! 

.  •   ■ 

MM!!!!!!!!'!!!!!!!!!!!!!!!!!!!'!!!!!!!!!!!! 

IH ill!    m 

Fig.  4.10:  Performance  of  coded  4-CPFSK  with  k=3 
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The  coding  starts  working  efficiently  after  6.7  dB  Eb/N0. 

At  this  signal-to-noise  ratio  almost  1  dB  gain  is  obtained  and 
at  7.6  dB  Eb/NQ,    the  coding  gain  reached  almost  2  dB. 

Signal-to-noise  ratios  for  8.8  dB  and  10  dB  were  also 
tested  and  no  errors  were  detected  in  fifty  thousand 
information  bits. 

The  MATLAB  source  code  for  this  system  is  listed  in 
Appendix  E . 

C.   CODED  4-CPFSK  WITH  CONSTRAINT  LENGTH  FOUR. 

The  final  coding  studied,  uses  four  shift  registers.  The 
encoder  is  illustrated  in  Fig.  4.11,  and  the  state  trellis  for 
this  system  is  illustrated  in  Fig.  4.12.  [Ref.  5] 

In  this  application,  the  demodulator  calculates  the 
euclidean  distances  of  the  received  signal  vector  from  eight 
different  vectors  presented  in  Table  4.4.  At  each  symbol 
duration  eight  distances  are  entered  to  the  Viterbi  decoder, 
and  the  true  path  is  chosen  according  to  minimum  distance. 

The  transition  matrix  of  this  system  is  illustrated  in 
Table  4.7  and  the  MATLAB  source  code  is  listed  in  Appendix  F. 

The  performance  of  this  coding  is  illustrated  in  Fig. 
4.13.  The  coding  starts  working  efficiently  at  5.9  dB  Eb/NQ. 

At  this  signal-to-noise  ratio  almost  1  dB  gain  is  obtained  and 
at  6.8  dB  Eb/NQI    the  coding  gain  reached  almost  2  dB. 
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Fig.    4.11:    Encoder. 


Most  significant  bit 


Least  significant  bit 


States 


Fig.  4.12:  Trellis  for  coded  4-CPFSK  with  k=4 . 
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TABLE  4.7:  VITERBI  TRANSITION  MATRIX. 


1  0 

4  3  7 

5  2  5 
7  3  7 

2  5 

3  8 


1102226225 


438313314 
52660260  1 


7  3 


1  4 


7  14  7  13 


1  3 


126202201 

3374  13414 


502501625626 


3  7  8  3 


Signal-to-noise  ratios  of  7.7  dB,  8.8  dB,  and  10  dB  are 
also  tested  and  no  errors  were  detected  for  fifty  thousand 
information  bits. 
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Fig.  4.13:  Performance  of  the  coded  4-CPFSK  with  k=4 
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V.  CONCLUSIONS  AND  FURTHER  RESEARCH 

This  thesis  mainly  presented  the  simulation  results  for 
rate  1/2  convolutional  quaternary  CPM  schemes  with  different 
constraint  lengths.  Using  different  codes,  the  minimum  merge 
length  increased,  and  larger  minimum  Euclidean  distances  were 
obtained.  Better  error  performances  were  obtained  as  a  result 
of  these  codes. 

This  research  can  be  extended  in  several  ways.  First,  rate 
1/2  codes  can  be  improved  by  using  different  constraint 
lengths,  and  shift  register  connections.  The  same  encoders  can 
also  be  tested  with  different  modulation  indexes. 

More  complicated  rate  1/2  codes  with  h=l/2,  can  be 
simulated  very  easily  by  making  small  changes  to  the  MATLAB 
source  codes  that  are  listed  in  Appendices  D  thru  F.  The 
signal  constellation  that  is  presented  in  Table  4.4  can  be 
used  for  every  new  code. 

Second,  rate  2/3  codes  for  different  constraint  lengths, 
shift  register  connections,  and  modulation  indexes  can  be 
studied,  and  best  systems  with  minimum  receiver  complexity, 
and  maximum  power -bandwidth  performances  can  be  found. 
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APPENDIX  A 


%This  program  creates  a  m  by  n  binary  input  matrix. 


function  [random_matrix] = random (m, n) ; 

a=rand (m; n) ; 

b=ones (m, n) *0  .  5  ; 

random  matrix=f loor (a+b) ; 


%    This  program  is  used  as  mapper.  It  converts  0  to  -1  and  1 
%  to  0. 

function  [mapper_output] =mapper (input) ; 
[m, a] =size ( input ) ; 
for  x=l:a 

y=ones (m, 1) *2^ (2-x) ; 

z= [z,y] ; 
end; 

mapper_output= input . *z; 
mapper_output=mapper_output ' ; 
mapper_output= (mapper_output) -1; 


%   This  program  is  used  as  CPFSK  modulator.  Inputs  to  this 
%  program  are  "T=symbol  duration",  "h=modulation  index",  and 
%  "modulator_input=mapper  output". 

function  [mod_output] =modulator (T, h,modulator_input) ; 

z=length (modulator_input ) ; 

teta0=0; 

for  sample=l:z 

I=modulator_input (sample) ; 

if  sample==l 

teta=0;  %Sets  the  initial  phase  to  zero. 

else 

tetaO=tetaO+modulator_input (sample- 1) ; 
teta=pi*h*tetaO; 

end 

t ime=0: 0.0005 :T- 0.00 05; 

time=time+ (sample- 1) *T; 

xx=pi*I*h* (time/T- (sample-1) ) +teta; 

mod_output (sample, : ) = (exp (i*xx) ) ; 
end 
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%    This  program  is  used  to  add  Gaussian  noise  to  modulator 
%  output.  Inputs  to  this  program  are  "m  as  modulator  output", 
%  "coeff  as  standard  deviation  of  the  noise",  and  "T  as  symbol 
%  duration".  The  constant  0.005  as  sampling  period. 

function  [awgn_noise] =awgn_ch (m, coeff , T) 

rand ( ' seed' , 99 ) ; 

rand ( ' normal ' ) ; 

awgn_noise=coeff * (rand (m, T/0 . 005 ) +j *rand (m, T/0 .005) ) ; 


%    This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
%  "mod_output  as  received  signal".  In  this  demodulator  hard 
%  decision  is  used. 

function  [demod_output] =demodulator (T, h,mod_output ) ; 
map= [ - 1  1]  ; 

[m,a] =size (mod_output) ; 
t ime=0: 0.0005 :T- 0.0005, • 
for  sample=l:m 
x=rem ( sample , 2 ) ; 
phase_coef f= [ (l-x) *pi] ; 
phase=exp (i*phase_coef f ) ; 
for  m_ary=l:2 

xx=exp (i*time*pi*h*map (m_ary) /T) ; 
match=xx*conj (mod_output (sample, : ) ' ) ; 
match=phase*match; 
match_out= [match_out  match] ; 
end 
[y, i] =max (real (match_out ) ) ; 
demod_output (sample) =map (ceil (i/1) ) ; 
match_out= [  ] ; 
end 
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function  [Nb,I,e]  =check(x,y) 

%  SYMBOL  ERROR  CHECK 

%  Paul  H.  Moose 

%  Naval  Postgraduate  School 

%  09-01-91 

%  This  m-file  locates  the  positions  in  vectors  x  and  y  that  do 

not 

%  agree.  It  returns  a  one  in  e  if  they  do  not  agree  and  a  zero 

if  they 

%  do  agree,  e  is  the  error  vector  if  x  and  y  are  binary.  I  is 

a  vector 

%  of  error  location  numbers.  Nb  is  the  sum  of  the  elements  of 

e . 

e= (x~=y) ; 

I=find(e) ; 

Nb=sum(e) ; 


%   This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system. 

clear 

T=0.015;  %  "T"  is  used  as  symbol  duration. 

h=l;      %  "h"  is  used  as  modulation  index. 

m=50000; 

n=l;      %  "m"  and  "n"  are  used  to  create  binary  matrix. 

coeff=0;  %  "coeff"  is  used  as  standard  deviation  of  the  noise. 

[input_matrix] =random(m;n) ; 

[mapper_output] =mapper (input_matrix) ; 

[mod_output] =modulato (T,h,mapper_output) ; 

[awgn_noise] =awgn_ch (m, coeff , T) ; 
mod_output=mod_output+awgn_noise; 

[demodulator_output] =demodula (T, h,mod_output) ; 
error=check (mapper_output , demodulator_output ) ; 
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APPENDIX  B 


%    This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
%  "mod_output  as  received  signal".  In  this  demodulator  hard 
%  decision  is  used. 

function  [demod_output] =demodulator (T, h,mod_output ) ; 
map= [-1  1] ; 

[m, a] =size (mod_output) ; 
t ime=0: 0.0005 :T; 
for  sample=l:m 

mod=mod_output (sample, : ) *exp ( j* (sample- 1) *pi/2) ;  %Of fsets 

%the  phase 
for  m_ary=l:2 

xx=exp  (i*time*pi*h*map  (m_ary)  /T)  ; 
match=xx*conj (mod_output (sample, : ) ' ) ; 
match_out= [match_out  match] ; 
end 
demod_output (sample, : ) = (real (match_out) ) ; 
match_out= [  ] ; 
end 

R=  ( j*demod_output  (  :  ,  1)  +demod_output  (:  ,2) )  /31; 
[demod_output]  =eucdis  (2  ,  R)  ;  %Calculates  the  euclidean  distance 

%of  the  received  signal 


function  D  = 

=  eucdis (q, R) 

% 

EUCLIDEAN  DISTANCE  METRICS 

% 

Paul  H.  Moose 

% 

Naval  Postgraduate  School 

% 

06-17-92 

%  This  M-file  finds  Euclidean  distance  of  elements 

%  in  vector  R  from  2^q  unit  amplitude  vectors 

%  equally  spaced  on  the  unit  circle.  It  stores  these  as  rows 

of  D. 

% 

N=2"q; 

L=length(R) ; 

index=l:N; 

dph=2*pi/N; 

M0=exp(j*dph.* ( index- 1)  )  ; 

for  1-1 :L 

D(l, : ) =abs(R(l) .*ones(M0) -M0) ; 

end 
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function  PHN  =  sof tv (k, K, Np, PH, T, D) 
%  Soft  Viterbi  Decoder 

%  Paul  H.  Moose 

%      Univ.  degli  Studi  di  Padova 
%  17-05-91 

%   This  M-file  decodes  k  bit  msgwords  from  2^n  real  metrics 

%  (These  may,  for  example,  represent  the   "distance"  of  the 

%  received  modulation  value  from  each  of  2^n  modulation 

values . ) 

%  The  state  transition  information  for  a  2AK  state  trellis  is 

in 

%  the  2AK  by  3*2^k  matrix  T.  Each  of  the  2Ak  entering  paths 

to 

%  each  state  has  its  source  state  (one  of  2^K) ,  path  msgword 

(one 

%  of  2Ak)  and  path  codeword  (one  of  2An)  listed  in  the  state 

row. 

%    The  path  histories  are  kept  in  matrix  PH  that  is  2^K  by 

3*Np. 

%  The  path  history  for  each  state  contains  source  state,  path 

%  weight  and  path  codeword  for  Np  previous  states. 

%   The  output  PHN  is  the  update  of  PH,  the  new  path  history. 

%  The  decoded  code  word  is  in  the  last  column  of  PHN.  (They 

should 

%  "merge" . 

%     The  past  histories  are  undated  on  the  basis  of  the 

"minimum 

%  metric".  You  can  change  this  to  the  "maximum  metric"  if 

desired  as 

%  indicated  in  the  comments  in  the  code. 

for  j=l:2AK 

X(j,2)=D(T(j  ,3)  ) +PH(T(j  ,1)  ,2)  ;  %path  weight 
X( j , 1) =T( j , 1) ;%path  source  state 

X ( j , 3 ) =T ( j , 2 ) ;  %path  code  word  T  ( j  ,  3  )  .  Chg  to  T  ( j  ,  2 )  for 
msgword. 

for  l=2:2^k 

wt  =  D(T(j,3*l))  +PH(T(j ,3*1-2) ,2) ; 

if  wt  <  X(j,2)  %The  <  selects  min  metric 

X ( j , 2 ) =wt ; 

X( j,l)=T(j, 3*1-2 ) ; 

X(j ,3)=T(j ,3*1-1) ;   %Chg   to   T(j,3*l)   for 
codeword. 

else 

end 

end 
%  We  now  need  to  append  old  paths  to  new  paths  to  get 
survivors . 

PHN(j,:)=  [X(j,:)  PH (X ( j , 1) , 1 : 3*Np- 3 ) ] ; 
end 
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%   This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  other  functions  are  same  as  in  Appendix  A. 


clear 

temp=[l  0  12  1  4;1  2  0  3];  %Is  used  as  Viterbi  path 

%matrix. 
T=0.015 
h=l/2; 
m=50000 
n=l; 
coeff=0 

[input_matrix] =random(m,n) ; 

[mapper_output] =mapper (input_matrix) ; 

[mod_output] =modulato (T, h,mapper_output ) ; 

[awgn_noise] =awgn_ch (m, coef f ,  T)  ; 
mod_output=mod_output+awgn_noise; 

[demodulator_output] =demodula (T, h,mod_output) ; 
viterbi_path_matrix=temp; 
TT=zeros (2, 60) ; 
for  x=l:m 

D=demodulator_output (x, : ) ; 

[TT] =softv(l, l;20,TT,viterbi_path_matrix,D) ; 
received_signal (x) =TT(1, 60) ; 
end 

input_matrix=input_matrix(19 :m) ; 
received_signal=received_signal (1 :m-19) ; 
error=check(received_signal, input_matrix) ; 


35 


APPENDIX  C 

%  This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 
%  "in  as  received  signal".  In  this  demodulator  hard  decision 
%  is  used. 

function  [demodulator_output] =demodula (T, h, in) 
N=T/.0005; 
[k, a] =size ( in) ; 
p  =  1 :  a ; 
y3— 1; 
for  m=3 :k 

phil=exp(j*pi* (m-l)/2) *exp(j* (pi/2) .* (p/N) ) ; 
phi2=exp(j*pi* (m-2) /2) *exp(j* (pi/2) .* (p/N) ) ; 
x=rem(m, 2 ) ; 
if  (x==0) 

y=sum(real (in(m-2, : ) ) .*abs (imag(phi2) ) )  + 

sum (real (in(m-l, : ) ) . *abs (imag (phil) ) ) ; 
if    (  (y<0    Sc  y3>0)     |     (y>0    &  y3<0)  ) 

d=l; 
else 

d=  - 1  ; 
end 
else 

y  =  sum(imag(in(m-2,  :  )  )   .*abs  (real  (  ph  i  2 )  )  )  + 

sum (imag (in(m-l, : ) ) . *abs (real (phil) ) ) ; 
if    (  (y>0    Sc  y3<0)     |     (y<0    &  y3>0)  ) 

d=  - 1  ; 
else 

d-1; 
end 
end 
y3=y; 

demodulator_output (m) =d; 
end 
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This  program  is  used  as  master  program.  It  calls  different 
functions  and  simulates  whole  communication  system.  All 
other  functions  are  same  as  in  Appendix  A. 


clear 
T=0.015 
h=l/2; 
m=50000 
n=l; 
coeff=0 
[input_matrix] =random(m,n) ; 
[mapper_output] =mapper (input_matrix) ; 
[mod_output] =modulato (T, h,mapper_output ) ; 
[awgn_noise] =awgn_ch (m, coef f , T) ; 
mod_output=mod_output+awgn_noise; 
[demodulator_output] =demodula (T, h,mod_output ) ; 
error=check (mapper_output (2 :m-2 ) , demodulator_output (4 :m] 
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APPENDIX  D 

%    This  program  simulates  the  encoder  that  is  given  in 

%  Fig.  4.4.  The  input  to  this  program  is  a  m  by  1  binary 

%  input  matrix.  At  the  end  it  gives  a  coded  m  by  2  binary 

%  output  .  First,  two  shift  registers  are  initialized  and  then 

%  input  are  shifted. 

function  [tcm_output] =tcm( input_matrix) ; 

[m, a] =size ( input_matrix) ; 

f f 1= [0; input_matrix( : , 1) ] ;  %  Intializes  the  first  flip-flop; 

ffl=ffl (l:m,l) ; 

ff2=  [0; ff 1]  ;  %  Intializes  the  second  flip-flop; 

ff2-ff2  (l:m,l)  ; 

tcm_output ( : , 1) =f f 1;        %  Most  significent  bit. 

tcm_output ( : , 2) =f f2 ;        %  Least  signif icicent  bit. 


%   This  program  is  used  as  mapper.  It  simulates  the 
%  Table  4.1. 

function  [mapper_output] =mapper (input ) ; 
[m, a] =size (input) ; 
for  x=l:a 

y=ones (m, 1) *2A (3-x) ; 

z= [z,y] ; 
end 

mapper_output=input . *z; 
mapper_output=mapper_output ' ; 
mapper_output  =  sum (mapper_output ) -  3 ; 
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%    This  program  is  used  as  CPFSK  demodulator.  Inputs  to  this 
%  program  are;  "T  as  symbol  duration",  "h  as  modulation  index" 

%  "mod_output  as  received  signal". 

function  [demod_output] =demodulator (T, h,mod_output ) ; 
mapl= [-3  -1  3] ; 
[m, a] =size (mod_output ) ; 
time=0 : 0.0005 :T; 
for  sample=l:m 

mod=exp ( j  * (sample- 1) *pi/2 ) *mod_output (sample ,  : )  ; 
for  m_ary=l:4 

xx=exp (i*time*pi*h*mapl (m_ary) /T) ; 
matchl=xx*conj  (mod' ) ; 
match_out= [match_out  matchl']; 
end 
demod_output (sample, : ) =distance (real (match_out) ) ; 
match_out= [] ; 
end 


%    This  program  calculates  the  euclidean  distances   of  the 
%  received  signal  from  signal  constellation. 

function  [dist] =distance (R) ; 

matrix=[0  -1  0  -31        %  This  matrix  represents  the  signal 

0  1  0  31  %constellation. 

1  0  31  0 
-1  0  -31  0 
0  -31  0  -1 
0  31  0  1 
31  0  1  0 

-  3 1  0  - 1  0  ]  ; 
d=ones (8,1)  ; 
dd=d*R; 

dd= (dd-matrix) . A2 ; 
dist=sum(dd' ) ; 
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%   This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 

%  other  functions  are  same  as  in  Appendix  A. 

clear 

temp=[l  01102213214 

125126237238]; 
T=0.015; 
h=l/2; 
m=50000; 
n-1; 
coef f =0; 

[input_matrix] =random(m,n) ; 

[tcm_output] =tcm(input_matrix) ; 

[mapper_output] =mapper (tcm_output ) ; 

[mod_output] =modulato (T, h,map) ; 

[awgn_noise] =awgn_ch (m, coef f , T) ; 
mod_output=mod_output+awgn_noise; 

[demodulator_output] =demodula (T, h,mod_output ) ; 
viterbi_path_matrix=temp; 
TT=zeros (2, 12  0) ; 
for  x=l:m 

D=demodulator_output (x, : ) ; 

[TT] =softv(2, 1, 40, TT, viterbi_path_matrix, D) ; 
received_signal (x) =TT(1, 120) ; 
end 

ww=mb (2, received_signal) ; 
ww=ww(l:2 :2*m) ; 
error=check(input_matrix(l:m-41) ' , ww(42 :m) ) ; 
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APPENDIX  E 


%   This  program  simulates  the  encoder  that  is  given  in 

%  Fig.  4..  The  input  to  this  program  is  a  m  by  1  binary 

%  input  matrix.  At  the  end  it  gives  a  coded  m  by  2  binary 

%  output  .  First,  the  three  shift  registers  are  initialized 

%  and  then  input  is  shifted. 

function  [tcm_output] =tcm(input_matrix) ; 

[m, a] =size ( input_matrix) ; 

f f 1= [0 ; input_matrix ( : , 1) ] ;  %  Initialize  the  first  flip-flop. 

ffl-ffl (l:m,  1) 

%  Initialize  the  second  flip-flop 


ff2=[0;ffl] ; 

ff2=ff2 (l:m,l) 

ff3=[0;ff2] ; 

ff3=ff3 (l:m;l) 

tcm_output ( : , 1) =abs (abs (f f l-f f2) -ff3)  ;  %  Most  significant  bit 

tcm_output ( : , 2) =f f2;  %  Least  significant  bit 


%  Initialize  the  third  flip-flop. 
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function  [PHN] =  viterbi (k, Np, PH, T, D) 

%       Soft  Viterbi  Decoder 

%         Paul  H.  Moose 

%      Mercury  Digital  Communications 

%  06-09-91 

%   This  M-file  decodes  k  bit  msgwords  from  2^n  real  metrics 

%  (These   may,  for  example,  represent  the   "distance"  of  the 

%  received  modulation  value  from  each  of  2^n  modulation 

values . ) 

%  The  state  transition  information  for  a  2AK  state  trellis  is 

in 

%  the  2^K  by  3*2^k  matrix  T.  Each  of  the  2Ak  entering  paths 

to 

%  each  state  has  its  source  state  (one  of  2AK) ,  path  msgword 

(one 

%  of  2^k)  and  path  codeword  (one  of  2An)  listed  in  the  state 

row. 

%    The  path  histories  are  kept  in  matrix  PH  that  is  2^K  by 

3*Np. 

%  The  path  history  for  each  state  contains  source  state,  path 

%  weight  and  path  codeword   for  Np  previous  states. 

%   The  output  PHN  is  the  update  of  PH,  the  new  path  history. 

%  The  decoded  msg  word  is  in  the  last  column  of  PHN.  (They 

should 

%  "merge". 

%     The  past  histories  are  undated  on  the  basis  of  the 

"minimum 

%  metric".  You  can  change  this  to  the  "maximum  metric"  if 

desired  as 

%  indicated  in  the  comments  in  the  code. 

a, 

o 

P=PH( : ,2) ' ; 
wt=T(: ,3:3:3*2"k) ' ; 
ux=T( : , l:3:3*2*k-2) ' ; 
aa=D(wt (:) ) ; 
bb=P(ux(:) ) ; 

wt ( : ) =aa+bb;  %This  contains  all  weights (columns)   for  each 
state (row) 

[a,b] =min(wt) ;  %  Use  max(wt)  here  for  maximum 
X(: ,2)=a' ; 

X( : , l)=diag(T( : ,3.*b-2) ) ; 

X( : ,3) =diag(T( : , 3 .*b-l) ) ;%path  msgword.   Chg   to   3 . *b   for 
codeword 

%We  now  need  to  append  old  paths  to  new  paths  to  get 
survivors . 
PHN= [X  PH(X( : ,1) ,l:3*Np-3) ] ; 
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%   This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  other  functions  are  same  as  in  Appendix  A. 

clear 

temp=[l  01102225226 

337338413414 

202201126125 

31431343843  7]; 
T=0.015; 
h=l/2; 
m=50000; 
n=l; 

coeff=0; 

[input_matrix] =random(m,n) ; 
[tcm_output] =tcm(input_matrix) ; 
[mapper_output] =mapper (tcm_output) ; 
[mod_output] =modulato (T, h,map) ; 
[awgn_noise] =awgn_ch (m, coef f , T) ; 
mod_output=mod_output+awgn_noise; 
[demodulator_output] =demodula (T, h;mod_output ) ; 
viterbi_path_matrix=temp; 
TT= zeros (4, 120) ; 
for  x=l:m 

D=demodulator_output (x, : ) ; 

[TT] =viterbi (2 , 40 , TT, viterbi_path_matrix, D) ; 
received_signal (x) =TT(1, 120) ; 
end 

ww=mb (2 , received_signal) ; 
ww=ww (1:2: 2*m) ; 
error=check( input (l:m-41) ' , ww(42 :m) ) ; 
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APPENDIX  F 


%    This  program  simulates  the  encoder  that  is  given  in 

%  Fig.  4.  The  input  to  this  program  is  a  m  by  1  binary 

%  input  matrix.  At  the  end  it  gives  a  coded  m  by  2  binary 

%  output  .  First,  the  four  shift  registers  are  initialized 

%  and  then  input  is  shifted. 

function  [tcm_output] =tcm(input_matrix) ; 
[m, a] =size ( input_matrix) ; 
ffl= [0;input_matrix( :,!)]; 


ffl=ffl(l:m 
ff2=[0;ffl] 
ff2=ff2 (l:m 
ff3=[0;ff2] 
ff3=ff3 (l:m 
ff4=  [0;ff3] 
ff4  =  ff4  (l:m 
tcm_output ( 
tcm_output ( 


1) 

1) 

1) 

1) 

,1) 

,2) 


=abs(ffl 
=  ff3: 


abs(ff2-ff4) ) 
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%    This  program  is  used  as  master  program.  It  calls  different 
%  functions  and  simulates  whole  communication  system.  All 
%  other  functions  are  same  as  in  Appendix  A. 

clear 

temp= [101102226225 

437438313314 

525526602601 

737738814815 

125126202201 

338337413414 

502501625626 

714713837838]; 
T=0.015; 
h=l/2; 
m=50000; 
n=l; 

coeff=0; 

[input_matrix] =random(m, n) ; 
[tcm_output] =tcm(input_matrix) ; 
[mapper_output] =mapper ( tcm_output) ; 
[mod_output] =modulato (T,h,map) ; 
[awgn_noise] =awgn_ch (m, coef f , T) ; 
mod_output=mod_output+awgn_noise; 
[demodulator_output] =demodula (T, h,mod_output) ; 
viterbi_path_matrix=temp; 
TT=zeros (8, 12  0) ; 
for  x=l:m 

D=demodulator_output (x, : ) ; 

[TT] =viterbi (2 , 40, TT, viterbi_path_matrix, D) ; 
received_signal (x) =TT(1, 120)  ; 
end 

ww=mb (2, received_signal) ; 
ww=ww(l :2 :m*2) ; 
error=check (input_matrix (1 :m-42) ' , ww(43 :m) ) ; 
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